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ABSTRACT 


Studies  of  high-pressure  ablation  mechanisms  have  been  hampered 
by  lack  of  suitable  ablation  models.  A  technique  has  been  developed 
to  contour  wedge  models  for  testing  in  the  Cornell  Aeronautical  Labs, 
wave  superheater  to  give  either  a  constant  pressure,  a  constant  heat 
transfer,  or  a  constant  pressure  gradient  up  to  the  limit  of  the 
facility.  Preliminary  testing  with  a  standard  carbon  phenolic  com¬ 
posite  shows  that  high  recession  rates  occur  even  in  the  absence  of 
a  significant  pressure  gradient. 
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NOMENCLATURE 


H  Enthalpy,  Btu/lb 

M  Mach  number 

p  Pressure,  atm 

o 

q  Cold  wall  heat  flux,  Btu/ft  -sec 

s  Distance  along  wetted  surface  from  holder  stagnation  point, 

in.  or  ft 

T  Temperature,  °R 

t  Time,  sec 

t^  Time  of  initial  exposure  of  the  model  to  the  flow  measured  from 

a  reference  time  zero  to  the  time  of  initial  flow  or  the  time 
the  teflon  cap  blew  off 

At  Time  a  model  is  exposed  to  the  flow,  sec 

U  Velocity,  ft/sec 

X  Distance  from  the  rotor,  in. 

x  ,  y  Orthogonal  contour  coordinates  with  origin  at  the  sample  leading 
8  edge,  in. 

a  Angle  of  attack 

3 

p  Density,  lb/ft 

Subscripts 

e  Local  conditions  just  outside  boundary  layer 

s  Gas  conditions  at  model  stagnation  point,  i.e.,  stagnation 

conditions 

00  Free  stream  gas  conditions 

Superscript 

o  Total  or  reservoir  conditions 
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I .  INTRODUCTION 


An  understanding  of  the  mechanisms  by  which  thermo chemical  processes  couple 

with  high-pressure  flow  fields  to  produce  excessive  ablation  rater  is  necessary 

before  reliable  analytical  procedures  can  be  developed  for  sizing  nose  tips  and 

flaps.  The  development  of  Improved  ablative  tip  materials  is  also  hampered  by 

our  lack  of  understanding  of  these  mechanisms.  Among  the  major  experimental 

difficulties  has  been  the  problem  of  characterizing  the  local  environment  both 

mechanically  and  thermally.  The  ablation  models  most  frequently  used  were  con- 

1  2 

leal  tips  of  various  nose  radii  and  half  angles,  45-deg  flat  wedges,  and 

3 

inclined  cylinders.  Both  of  the  latter  employed  water-cooled  leading  edges. 

In  all  these  cases,  and  especially  in  the  tip  models,  large  and  constantly 
changing  pressure  gradients  along  the  model  axis  made  interoretation  of  the 
results  difficult.  This  was  seriously  aggravated  by  the  fact  that  the  models 
changed  shape  during  the  test  and  the  environment  became  even  more  difficult 
to  define. 

Nevertheless,  there  was  a  high-pressure  gradient,  and,  for  tip  models,  the 

1  2 

rate  of  removal  could  be  correlated  with  stagnation  pressure.  *  As  a  conse¬ 
quence,  several  plausible  theories  were  evolved  in  which  pressure  gradient  shear 

2  A 

was  the  driving  mechanism.  ’ 

In  addition  to  the  macroscopic  pressure  gradient,  it  was  postulated  that 
surface  roughness  produced  by  inhomogeneous  ablation  (due  to  local  difference 
in  resin  content,  fiber  content  and  density)  was  coupling  with  the  boundary 
layer  to  produce  large  local  forces  and  heat  transfer  rates,  which  could  be 
responsible  for  the  massive  ablation  rates  observed. 

It  became  apparent  that,  in  order  to  get  a  quantitative  estimate  of  how 
the  various  processes  and  material  variables  were  interacting,  it  would  be 
necessary  to  develop  Improved  ablation  models.  The  ideal  model  would  have  a 
relatively  large  area  to  minimize  or  control  the  effect  on  resin-reinforcement 
geometries  and  would  be  shaped  to  control  pressure,  pressure  gradient,  and  heat 
transfer  so  that  each  of  these  variables  could  be  held  constant  over  the  entire 
model  surface.  It  was  with  this  objective,  therefore,  that  the  following  work 
was  undertaken. 
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Table  1.  Test  Conditions 


II.  TEST  FACILITY 


The  Cornell  Aeronautical  Labs.  (CAL)  wave  superheater  test  facility  was 
chosen  because  of  the  high  pressures  and  relatively  large  flow  field  that  are 
available.  This  hot  gas  generator  is  a  rotating  multiple-shock-tube  device 
capable  of  delivering  a  5-lb/sec  flow  of  uncontaminated  test  gas  with  effective 
reservoir  conditions  up  to  2800-Btu/lb  enthalpy  and  130-atm  pressure  for  dura¬ 
tions  up  to  15  sec.  Energy  is  transferred  from  the  high-pressure,  preheated 
test  gas  by  a  shock  wave  mechanism  in  the  288  shock  tubes  mounted  on  the 
periphery  of  a  rotor.  The  resulting  steady  flow  from  the  rotor  is  supersonic, 

with  height  and  width  dimensions  of  approximately  1.5  and  0.5  in.  The  high 

2 

power  density  (up  to  27,000  hp/in.  )  of  the  freely  expanding  jet  is  ideally 
suited  for  high  heat  transfer,  pressure,  and  shear  stress  testing.  The  flow 
from  the  free  jet  is  inherently  unsymmetrical ;  however,  this  effect  was  mini¬ 
mized  in  the  model  design.  For  this  study  we  selected  a  stagnation  pressure 

of  80  atm  so  that  the  ablation  data  could  be  compared  with  those  of  earlier 
2 

experiments.  Details  of  the  hot  gas  environment  are  given  in  Table  1. 
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Figure  1.  Water-Cooled 
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III.  MODEL  DEVELOPMENT 


Following  a  recommendation  by  Tate,^  we  decided  to  start  with  the  standard 
45-deg  wedge  model  with  a  water-cooled  leading  edge  (Fig.  1)  and  attempt  to 
contour  the  model  surface  to  obtain  constant  pressure  gradients.  This  wedge 
has  a  number  of  advantages:  it  has  a  large  surface  area,  1.5  in.  wide  x  2.0  in. 
long;  mounting  the  specimen  so  that  the  1.5-in.  dimension  is  vertical  nullifies 
the  effect  of  assymetries  in  the  flow  field  and  enhances  the  conditions  on  the 
model;  the  wedge  has  a  fixed  leading  edge. 

During  wedge  tests,  the  flow  is  stagnated  on  the  water-cooled  leading  edge 
of  the  specimen  holder.  The  gas  flows  from  this  region  over  the  specimen. 

If  the  pressure  is  constant  or  decreasing,  the  shock  systems  are  weak  and  the 
flow  over  the  specimen  is  nearly  isentropic.  In  this  case,  the  holder  stagna¬ 
tion  pressure  is  also  the  total  pressure  of  the  flow  over  the  specimen,  and  the 
ratio  of  measured  local  static  pressure  to  holder  stagnation  pressure  yields 
the  local  Mach  number  through  the  isentropic  flow  relationship.  The  total 
enthalpy  of  the  flow  outside  the  boundary  layer  is  constant,  since  heat  is  lost 
only  in  the  thin  boundary  layer.  From  the  local  pressure  distributions  presented 
in  the  following  sections,  the  rest  of  the  flow  conditions  over  the  specimen  can 
be  calculated.  The  total  conditions  of  the  flow  from  the  rotor,  the  holder  stag¬ 
nation  conditions,  and  static  conditions  0.5  in.  from  the  specimen  leading  edge 
are  presented  in  Table  1. 

Our  first  attempt  was  to  develop  a  constant-pressure  model.  In  order  to 
minimize  the  pressure  gradient,  we  used  two  methods  to  calculate  a  curvature 
for  the  specimen.  One  method  is  a  weak-shock-wave  solution  technique  based  on 
the  Prandtl-Meyer  turning  angle;  it  is  described  in  the  Appendix.  The  other  is 
a  strong-shock-wave  solution,  which  incorporated  modified  Newtonian  theory. 

Both  methods  make  maximum  use  of  flow  field  calibration  data  and  the  45-deg 
flat  wedge  pressure  distribution  obtained  in  earlier  programs. 

We  obtained  a  reasonably  high  pressure  over  a  long  test  surface  by  setting 
the  initial  angle  of  the  specimen  at  35  deg,  with  the  model  holder  stagnation 
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point  0.5  in.  from  the  rotor.  The  constant-pressure  contours  calculated  by 
the  Prandtl-Meyer  and  modified  Newtonian  methods  were  significantly  different, 
and  a  compromise  contour,  halfway  between  the  two  calculated  contours,  was 
selected  for  the  first  calibration  model  (Fig.  2).  Solid  copper  calibrating 


COORDINATE  PARALLEL  WITH  SPECIMEN  LEADING  EDGE  (xs),  in. 

Figure  2.  Comparison  of  the  Model  Contours 

models  were  instrumented  with  both  calorimeters  and  pressure  ports,  as  shown 
in  Fig.  3  and  Fig.  A,  which  is  a  photograph  of  a  mounted  model.  A  cap  con¬ 
sisting  of  two  1/8-in. -thick  strips  of  Teflon  was  wrapped  around  the  holder 
and  model  and  fastened  to  the  top  and  bottom.  This  cap  was  used  to  delay 
model  heating  until  the  flow  had  become  established.  The  pressure  distribution 
obtained  from  this  model  is  shown  in  Fig.  5.  The  centerline  pressures  seem  to 
indicate  an  overexpansion  near  the  leading  edge  and  a  strong  recompression 
toward  the  aft  portion  of  the  model.  The  initial  overexpansion  was  not  well 
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Figure  A.  Typical  Calibration  Model 

defined,  because  its  presence  was  indicated  by  only  one  pressure,  which  may 
have  been  influenced  by  the  presence  of  the  calorimeter  immediately  ahead  of 
it.  However,  the  well  defined  recompression  indicated  that  the  contour  was 
too  steep. 

The  shallower  contour  calculated  by  the  Prandtl-Meyer  method  was  selected 
for  the  second  attempt  at  a  constant-pressure,  model.  Tolerances  on  the  second 
model  were  held  to  limits  used  for  hypersonic  throats,  and  test  surface  coor¬ 
dinates  were  within  0.005  in.  of  the  designated  coordinates.  In  addition,  the 
calorimeters  and  pressure  taps  were  moved  to  separate  rows  to  reduce  any  pos¬ 
sible  interactions  (Fig.  2).  This  model  generated  a  nearly  constant  pressure 
distribution,  as  shown  by  the  dashed  curve  in  Fig.  6,  although  the  pressure 
level  was  about  25  percent  higher  than  expected.  Motion  picture  examination 
revealed  that  the  jet  from  the  rotor  caused  the  model  to  move  about  0.A  in. 
and  rotate  about  5  deg.  The  rotation  increased  the  angle  of  attack  of  the 
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PRESSURE  RATIO,  pe/ps 


DISTANCE  FROM  SPECIMEN  LEADING  EDGE  ,  in. 


Figure  5.  Pressure  Distribution  from  the  Compromise  Contour 

surface,  thus  increasing  the  pressure.  Also,  the  lateral  movement  caused  the 
jet  to  impinge  on  the  inclined  surface  without  stagnating  at  the  geometrical 
stagnation  point  of  the  holder,  which  may  have  contributed  to  the  high  pres¬ 
sure.  Two  additional  models  generated  pressures  of  the  same  level  as  those 
predicted  by  the  calculations.  These  measured  distributions  are  presented  in 
Fig.  6. 
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DISTANCE  FROM  SPECIMEN  LEADING  EDGE,  in. 


Figure  6.  Pressure  Distributions  on  the  Prandtl-Meyer 
Constant-Pressure  Contour 


The  Prandtl-Meyer  method  was  then  used  to  generate  a  contour  to  create  a 
pressure  gradient  of  -10  atm/in.  The  leading  edge  angle  was  increased  from 
35  to  40  deg  to  bring  the  average  surface  pressure  close  to  the  level  generated 
on  the  constant-pressure  surface.  The  pressure  distribution  was  not  quite 
smooth,  but  the  average  gradient  over  the  sample  length  was  very  close  to 
-30  atm/in.,  as  shown  in  Fig.  6. 

Pressure  distributions  along  the  flat  45-deg  wedge,  which  were  obtained 
2 

in  previous  studies,  are  also  shown  in  Fig.  6.  All  pressures  shown  in  this 
figure  have  been  normalized  with  the  wedge  holder  stagnation  pressure,  which 
is  the  effective  total  pressure  for  the  flow  adjacent  to  the  sample.  The 


J 
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normalized  pressure  allows  direct  calculation  of  the  local  Mach  number  and 
the  other  local  flow  parameters. 

Heat  flux  data  for  the  two  contoured  and  flat  configurations  are  presented 

in  Fig.  7.  The  data  for  the  constant  pressure  contour  represent  an  average  of 

data  obtained  during  the  three  runs  when  the  pressure  gradient  was  close  to 

zero,  and  the  data  for  the  flat  wedge  represent  an  accumulation  of  previous 
2 

data.  Data  for  the  three  configurations  are  reasonably  close  to  the  predic¬ 
tion  of  turbulent  boundary  layer  theory  for  a  flat  plate,  with  the  distance 
parameter  measured  from  the  water-cooled  holder  stagnation  point.  Aerodynamic 
shear  stress,  calculated  from  the  Reynolds  analogy,  is  also  shown  in  Fig.  7. 
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DISTANCE  FROM  HOLDER  STAGNATION  POINT  s,  in. 


Figure  7.  Heat  Flux  and  Shear 
Stress  Distributions 


IV.  ABLATION  TESTING 


In  order  to  test  the  usefulness  of  these  new  contoured  models,  we  planned 
a  series  of  ablation  runs  using  a  REST  6300  HP  composite.  This  material,  devel¬ 
oped  under  a  SAMSO  contract  for  the  REST  program,  has  been  extensively  charac¬ 
terized  and  tested.  It  has  been  used  as  a  standard  under  the  RESEP  program 
for  the  past  few  years.  The  specimens  prepared  for  this  program  contained 
44  percent  Monsanto  SC1008  phenolic  and  56  percent  3M  Pluton  B-l  high  purity 
fabric  and  were  substantially  within  the  specification  for  this  material  devel¬ 
oped  under  the  REST  program.  All  specimens  were  cut  from  a  flat  laminate  4  in. 
thick  and  approximately  1  ft  square.  The  standard  45-deg  wedges  were  cut  to  a 
fabric  orientation  angle  of  30  deg  to  the  surface  away  from  the  direction  of 
flow.  Since  both  contoured  models  (0  and  10  atm/in.  pressure  gradients)  could 
not  be  cut  to  a  constant  fabric  orientation,  an  average  orientation  of  about 
30  to  35  deg  at  the  midpoint  (fore  to  aft)  was  selected.  Copper  adaptor  blocks 
were  designed  as  shown  in  Fig.  8  to  accommodate  the  contoured  specimens  to  the 
45-deg  wedge  model  holder. 

Environmental  conditions  and  specimen  performance  data  are  listed  in 
Table  2.  The  most  striking  data  are  those  for  total  recession.  For  all  three 
configurations  ac  the  point  of  maximum  recession,  the  amount  of  erosion  is  the 
sam'-.  The  constant-pressure,  constant-heat-flux  model  eroded  uniformly  over 
its  entire  surface,  as  was  expected.  A  comparison  of  the  Initial  and  final 
contours  of  this  model  with  those  of  the  45-deg  wedge  is  shown  in  Fig.  9.  The 
very  high  recession  rate  in  the  absence  of  any  significant  pressure  gradient 
shear  is  most  surprising.  Inspection  of  film  taken  during  the  runs,  at  600 
frames/sec,  shows  very  definitely  that  char  particles,  some  of  them  quite  large, 
are  being  removed  by  what  appears  to  be  some  mechanical  process.  It  is  inter¬ 
esting  to  note  that  the  45-deg  flat  wedge  model,  as  it  ablates,  tends  to  be 
the  same  shape  as  the  constant-pressure  contoured  models.  Under  the  REST  pro¬ 
gram  (Ref.  7)  copper  calibration  models  furnished  with  pressure  ports  were 
made  to  mockup  the  shapes  produced  as  the  model  ablated.  At  the  point  of  maxi¬ 
mum  recession,  the  local  pressure  rose  rapidly  to  the  same  local  pressure  as 
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FLAT  SPECIMEN 


Figure  8.  Configurations  of  Wedge  Models 
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Table  2.  Typical  Test  Data  for  Standard  Carbon  Phenolic 


Model  No. 

1 

2 

3 

Configuration 

Flat  Wedge 

dp/ds  ■  0  atm/in. 

dp/ds  -  -10  atm/in. 

a,  deg 

45 

35 

40 

X  Station,  in. 

0.5 

0.5 

0.5 

t1>  sec 

0.73 

1.02 

0.83 

t  (Run  Time) ,  sec 

1.77 

1.58 

1.62 

Specimen  weight,  g 

Before 

42.31 

28.30 

30.50 

After 

26.80 

16.85 

17.50 

Max.  Recession,  in. 

0.2545 

0.2531 

0.2644 

Location  of  Max.  Recession 

x  ,  in. 
s 

0.6 

Uniform 

0.8 

Avg.  Surface  Temp.,  °R 

x  ■  0.5  in. 
s 

5800 

5800 

5800 

x  ■  1.0  in. 
s 


5500 


5650 


5650 


Figure  Comparison  of  Contoured  Ablation 
Specimen  and  Flat  Specimen 

that  which  occurs  over  most  of  our  constant  pressure  model  ('v  30  atm).  The 
recession  rates  of  both  the  0-atm/in.  model  and  the  20-atm/in.  model  were 
constant  during  the  entire  run  and  were  equal  where  the  local  pressures  were 
equal.  These  data  indicate  that  the  recession  correlates  more  with  local 
pressure  than  with  pressure  gradient. 
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V.  CONCLUSIONS 


Ablation  models  with  relatively  large  surface  areas  can  be  designed,  with 
the  use  of  the  Prandtl-Meyer  turning  angle  technique,  to  give  constant  pressure 
gradients  in  the  CAL  wave  superheater.  For  the  case  where  the  pressure  gra¬ 
dient  is  zero,  the  heat  flux  is  also  a  constant.  In  this  case,  the  aerodynamic 
shear  stress  as  calculated  from  the  Reynolds  analogy  is  also  constant. 

For  the  limited  amount  of  ablation  testing  performed  with  a  standard 
R6300  HP  carbon  phenolic,  high  ablation  rates  with  apparent  mechanical  removal 
occur,  even  on  a  zero-pressure-gradient  model.  A  comparison  of  this  specimen 
with  a  45-deg,  20-atm/in. -pressure-gradient  model  at  the  point  of  maximum 
recession  appears  to  correlate  recession  rate  with  local  pressure  rather  than 
with  pressure  gradient. 
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APPENDIX.  DESIGN  PROCEDURE  FOR  A  CONTOURED 
WEDGE  WITH  AN  ARBITRARY  PRESSURE  GRADIENT 


There  are  three  basic  steps  in  the  procedure  for  generating  a  wedge  con¬ 
tour  for  use  in  the  wave  superheater  free  jet.  These  are  problem  specification, 
model  orientation,  and  the  actual  calculation  of  the  contour.  The  first  of 
these,  the  problem  specification,  includes  such  items  as  pressure  level,  pres¬ 
sure  gradient,  type  of  material  to  be  tested,  and  model  size.  Several  param¬ 
eters  should  be  considered  before  the  model  is  designed.  Obviously,  the  highest 
model  pressures  can  be  obtained  if  the  wedge  configuration  is  positioned  close 
to  the  rotor.  However,  in  many  cases  the  same  surface  pressure  can  be  obtained 
at  different  test  locations.  The  stagnation  pressure  on  the  wedge  holder 
leading  edge  is  the  effective  total  pressure  of  the  flow  over  the  wedge.  As 
the  holder  is  moved  farther  from  the  rotor,  the  stagnation  pressure  decreases, 
and  the  holder  incidence  angle  must  be  increased  to  maintain  the  same  static 
pressure  over  the  surface.  The  increased  incidence  angle  results  in  reduced 
Mach  number,  heat  flux,  and  shear  stress.  Therefore,  specification  of  surface 
pressure  alone  Is  not  sufficient  to  define  the  contour. 

Once  the  problem  has  been  properly  specified,  the  next  task  is  to  pick  a 
rotor  station  and  determine  the  initial  angle  of  the  sample  necessary  to  gen¬ 
erate  the  correct  pressure  at  the  leading  edge.  This  angle  is  determined  by 
generating  an  expression  for  the  pressure  distribution  over  a  flat  wedge  with 
the  same  angular  orientation.  Since  calibration  data  are  limited  to  45-deg 
wedges  at  a  few  stations  from  the  rotor,  this  distribution  is  not  known  a 
priori .  It  is  assumed  that  the  static  pressure  distribution  in  the  free  stream 
is  representative  of  the  distribution  over  a  0-deg  wedge  configuration.  This 
pressure  distribution  is  written  in  the  form 


where  a  and  b  are  constants,  p  is  the  stagnation  pressure  on  the  wedge  holder, 

s 

and  x  is  the  distance  from  the  rotor.  The  measured  or  interpolated  pressure 

distribution  over  a  45-deg  wedge  is  converted  into  a  similar  expression.  It 

should  be  reiterated  that  pg  is  also  the  total  pressure  of  the  flow  over  the 

wedge  surface  and  that  Pe/Ps  specifies  the  local  Mach  number.  For  the  cases 

considered  to  date,  the  exponent  b  is  nearly  the  same  for  both  the  O-deg  and 

the  45-deg  wedges.  The  desired  pressure  at  the  wedge  leading  edge  has  been 

specified,  p  is  determined  by  specifying  the  distance  from  the  rotor,  and 
s 

the  constant  a  can  be  determined.  The  variation  of  a  is  assumed  to  be  linear 
with  angle  from  0  to  45  deg,  and  the  initial  angle  of  the  specimen  is  obtained 
through  linear  Interpolation.  When  the  angle  is  known,  the  constant  b  can  be 
corrected  for  the  small  difference  between  0  and  45-deg  wedges,  if  necessary, 

and  the  calculation  repeated.  The  initial  angle  should  be  checked  for  consis¬ 

tency  with  the  problem  specifications.  At  this  point,  the  flow  conditions  at 
the  sample  leading  edge  and  the  distributions  on  a  flat  wedge  with  the  same 
leading  edge  angle  are  known. 

The  next  step  involves  the  actual  mechanics  of  generating  the  contour 
coordinates.  Simply  stated,  the  Mach  number  distribution  on  the  flat  wedge 
and  the  required  Mach  number  distribution  according  to  the  specified  pressure 
gradient  are  computed.  Then  the  Prandtl-Meyer  angles  for  both  distributions 
are  generated.  The  local  angles  of  the  flat  wedge  are  modified  by  the  dif¬ 
ference  between  the  specified  and  the  flat  wedge  Prandtl-Meyer  angles,  and  the 
contour  coordinates  are  obtained  by  a  finite-difference-type  procedure.  The 
step-by-step  calculation  procedure  is  as  follows 

1.  Select  a  nominal  value  for  Ax 

2.  Compute  x  for  each  value  of  Ax. 

3.  Compute  the  flat  wedge  values  for  pe/ps  =  ae  ^X,  Me  *  f(pe/ps), 
ve  *  f(Me),  and  Axs  =  Ax/cos  0O  for  each  value  of  x.  v  is  the 
Prandtl-Meyer  function,  and  a  calculation  method  for  real  gas  is 
given  by  Owczarek.l 


^Owczarek,  J.  A.,  Fundamentals  of  Gas  Dynamics ,  International  Textbook 
Company,  Scranton,  Pa.,  1964. 
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4.  Compute  Ays  Axa  tan  (0j_i  -  0O)  where  j  is  the  calculation 
index;  i.e.,  J  ■  0  at  the  sample  leading  edge  and  j  ■  1  at 
the  end  of  the  first  Ax8  Increment. 


5.  Compute  As,  ■  /  Ax  +  y 

J  8  H 

6.  Compute  the  specified  pressure: 


(!Wc/pS,j 


(p  /P  )<  n  +  ^  J  As. 
e  rs  j-0  ds  J 


7.  Compute  specified  Mach  no.  using  real  gas  functions: 


M  -  f(p  /p  ), 

spec  j  rspec  rs  j 


8.  Compute  the  specified  Prandtl-Meyer  function: 


v  . 
spec  j 


f ( M  _  J 

spec  j 


) 


9  .  Compute  the  new  wedge  angle : 


0o  + 


(v 


ej 


V  .  ) 

spec  j 


10.  A  slight  error  is  incorporated  in  Step  4  when  8^_p  the  angle 

of  the  previous  element,  is  used  as  the  angle  or  the  jth  element. 
This  error  is  usually  small  if  Ax  is  small,  but  it  can  be 
eliminated  for  all  practical  purposes  by  use  of  an  average  angle: 


Ays  -  Axg  tan[ (6j  +  0j _x> /2  -  eQ] 


Repeat  Steps  3  through  9  to  obtain  the  more  accurate  value  for 
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11.  Update  the  j  Index  and  start  from  Step  4  again.  Follow  this 
procedure  until  the  computed  value  of  As  ■  EAsj  exceeds  the 
specified  model  length  or  the  value  of  0j  exceeds  90  deg. 

12.  The  final  coordinates  of  xg  and  ys  should  be  smoothed,  since 
the  contour  is  as  sensitive  to  bumps  as  is  the  throat  region 
of  a  contoured  hypersonic  nozzle. 
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